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Abstract 

^ ■ We study the accretion of cold and hot dark matter onto a cosmic 

string wake using a high-resolution numerical simulation. We verify 
previous analytical calculations predicting the radius of bound mat- 
ter around wakes and inflow velocities of the dark matter, as well as 
assumptions about the self-similarity of the solution. In cold wakes, 
we show that self-similarity is approached quickly after a 'binding' 
transition. Hot wakes approach self-similarity rapidly once the free 
streaming 'pressure' falls below a critical value and accretion begins. 
We also analytically calculate the size of the overdensity in wakes with 
cold dark matter and compare the results to simulations. We remark 
that the results derived may be used in testing gravitational numerical 
codes in the non- linear regime. 



1 Introduction 

Large-scale structure formation in the cosmic string model is not yet well 
understood, particularly in the non-linear regime. Until recently, only linear 
analytical treatments have been available for approximating wake character- 
istics. Rees 1=1 has analytically calculated some of the non-linear effects of 
baryons in wakes, particularly at the time of formation. 

There have also been some recent numerical attempts at understand- 
ing the matter distribution better. Shellard and Avelino o have developed 
Zel'dovich techniques including the source terms from cosmic strings, but 
this work is only valid on large scales where the perturbations are still linear. 
Sornborger, et. al. o have simulated non-linear evolution of baryons and 
dark matter in cosmic string wakes, but their simulations had small grids, 
and therefore could not resolve non-linear features in the matter distribution 
well. Their simulations also have the drawback that boundary conditions 
come into play typically by a redshift of 8, and therefore assumptions have 
to be made about the scaling of the matter distribution (i.e. that it is self- 
similar). 

Therefore, what is needed is a confirmation of both linear and non-linear 
results and evolution of the wakes until the present in order to verify that the 
assumptions about the growth of the matter distribution in the wake are cor- 
rect. We would like more accurate analytical expressions for the size of wake 
overdensities to compare with observations. Also, since cosmic strings with 
cold dark matter generate too much small scale structure Q , an investigation 
of the effects of the accretion of hot dark matter would be desirable. 

In this paper, we have done high-resolution one-dimensional numerical 
simulations yielding precise comparisons with the Zel'dovich approximation 
in the linear regime. In addition, we have extended the Zel'dovich approxi- 
mation, in a slightly modified guise, to make predictions for the width of the 
overdensity in the wake into the non-linear regime. We have also extended 
our simulation to include the effects of hot dark matter. 

The cosmic string model of large-scale structure is motivated by the fact 
that supermassive strings are a common product of a GUT scale phase tran- 
sition. A string network is formed at the phase transition. The energy carried 
by the defects curves spacetime, and thus seeds structure formation. Simple 
arguments "t$, verified by numerical simulations, show that the string net- 
work rapidly evolves to a scaling solution. That is, the number of long strings 



per horizon volume is statistically constant (i.e. when we average over many 
volumes) . 

Numerical simulations have shown that cosmic string loops have a typical 
radius at the limit of resolution of cosmic string network simulations R ~ 
10 _2 i/ _1 , whereas there are ~ 10 — 30 horizon sized long strings per particle 
horizon volume in the matter era O'lil'O. This indicates that the primary 
structure formation mechanism in the cosmic string model is in the form of 
sheetlike wakes formed in the path of long strings. 

The matter perturbation caused by a moving long cosmic string comes 



from two conceptually distinct components of the gravitational field eU'liU'lii 
a velocity kick imparted to particles toward the plane traversed by the cos- 
mic string due to converging geodesies in the conical spacetime around the 
string, and a Newtonian field due to the effective mass of small-scale struc- 
tures such as waves and kinks moving on the string. In linear theory, the 
comoving velocity impulse toward the wake imparted to a particle after a 
string has passed by is 

2itGU efJ -T) A „ 
Vkick = 4:7rGfi eff v s -f s 

where G is Newton's constant, fi e ff is the effective mass-per-length of the 
string, v s is the velocity of the string, 7 S is the relativistic gamma factor 
7 S = (1 — f 2 ) -1 / 2 and T is the tension on the string which obeys an effective 
equation of state fi e ffT = //o- For strings formed at a GUT phase tran- 
sition, G/J.Q ~ 10~ 6 . In the expression for Vkick, the first term comes from 
the Newtonian field and the second from the deficit angle of the spacetime. 
Thus, we can see that for fast moving strings (v s ~ 1) the impulse from the 
Newtonian field is negligible compared to that from the deficit angle in the 
conical spacetime around the string. 

To study matter accretion in cosmic string wakes, planar symmetry is 
often assumed as a first order approximation for ease of calculation. That is, 
we assume the wake to be formed by a cosmic string in the regime where the 
string velocity is relativistic (a good approximation) and the string straight 
(a good approximation if we neglect small scale structure on the strings, 
which is at present not well understood O'lil'O, but can be hoped to be 
of the order of the cosmic string loop scale < 10~ 2 H). Thus, we neglect 
the Newtonian field of the string, and we assume that the string has passed 
through the volume of interest without deviating from being straight, such 



that the inflow is effectively one-dimensional. In this scenario, two matter 
streams are formed about the axis where the cosmic string has passed, each 
with identical velocities toward the axis. 

Analytical calculations using the Zel'dovich approximation tLrMLif and 



assuming planar symmetry have shown that the wake width should grow 
as £ 2 / 3 in the matter dominated era. In these calculations, the comoving 
wake width was defined as the physical radius h from the wake at which the 
particle turns around (h = 0). For hot dark matter, the width of the wake 
has been calculated including the effects of the larger phase space volume 
occupied by dark matter. These calculations predict that hot dark matter 
wakes laid down at a t eq will be of a size the same order as that of cold dark 
matter wakes laid down at the same time. 

Since the perturbation overdensity Sp/p = 1 in the wake immediately 
following passage of the string, the solution is non-linear ah initio inside the 
overdensity. For times t 3> ti, where t, is the formation time of the wake, it 
can be shown that matter accretion has no preferred length or time scales, 
thus the matter distribution in the wake approaches self-similarity 

In this paper, we study planar symmetric wakes numerically and derive 
and compare where we can with analytical results. We have developed a 
high-resolution one-dimensional particle-mesh code to study hot and cold 
dark matter accretion in cosmic string wakes in the non-linear regime. The 
calculations are done assuming a matter dominated background. 

In the next section, we analytically calculate the size of the overdensity in 
the wake. Then, in section 3, we present details of the initial and boundary 
conditions and methods used for the simulations. In section 4, we give re- 
sults of the numerical simulation and compare these with our own and other 
previous analytical results. Finally, we discuss our results and present our 
conclusions. 



2 The Width of the Wake Overdensity 

The equations in Lagrangian coordinates for a dark matter particle interact- 
ing with a gravitational field in an expanding universe in one-dimension are 



(see for instance [[U|] ) : 



v + 2-v = — M (1) 

a a 6 



r = v (2) 

d 2 <j) = 4nG{p-p). (3) 

We are working in comoving coordinates where v = (l/a)v p — (a/a)x p is 
the comoving particle velocity (subscript p indicates physical quantities), 
a = a(t) is the scale factor, cf> = acf) p + (a 2 /2)aXp is the comoving gravita- 
tional potential, r = r p /a is the comoving particle position, p = a 3 p p is the 
comoving mass density, and p is the average comoving background density. 
An overdot indicates derivative with respect to time, d denotes a spatial 
derivative and we normalize the scale factor to unity at the present time. 

Using the Zel'dovich approximation, it has been calculated that the co- 
moving infall velocity of particles outside the wake is (see, for instance 113) 

tw™ =§«,(£) 3 +H£) (4) 

where subscript i indicates the initial value of the variable (i.e. wake forma- 
tion). 

As mentioned in the introduction, analytical calculations using the Zel'- 
dovich approximation have estimated the width of the wake by calculating 
the radius at which the physical velocity of a particle goes to zero with respect 



to the plane of the wake iirl. However, the overdensity is bounded by the 
radius at which particles turn around in comoving coordinates, thus forming 
the first caustic in the density distribution. We refer to this radius as the 
'second turnaround radius'. We can calculate this radius as follows: 

We study a particle lying initially at the symmetry axis just as it is kicked 
by the passage of a string. We follow the particle's trajectory after the kick 
occurs. A particle at the axis when the kick occurs will travel away from the 
symmetry axis and be at the front edge of the overdensity. Thus the entire 
overdensity is behind it and the density can be given as a time dependent 
part times a delta-function distribution: 

p = p + p(t)S(x) (5) 

We can then reduce equations [1], H and [| to 

(a 2 v) = T-2nGp(t) (6) 

a 



= 2vrGp(t)|x| 



(7) 
(8) 



where — indicates a particle moving away from the wake and + is a particle 
moving towards the wake. 

We know the rate at which the particles are falling into the wake from 
equation |. Thus, the overdensity is growing as pit) = 2pbX str . eam (t), where 
Xstream is the distance particles have fallen towards the wake, given by the 
stream velocity integrated from the initial time to time t. That is 



/ N 6 

p[t) = -Vip b t 



(9) 



Using eq. |, and p b = l/6nGtl and v kickp = v kick = Vi(U/t ) 2/3 , where t is 
the present time, and ti is the time of wake formation, we find the comoving 
particle velocity 



v(t) = v kick 
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which is consistent with that found from the Zel'dovich approximation; and 
the particle trajectory 



x(t) = v kick 
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These results are exact until the particle's trajectory crosses with that of 
another particle. 

Second turnaround occurs at t — t turn , given by v{t turn ) = 0. Using eq. 
|17j| we find t tur n = bti, where b = 4.452. This gives a second turnaround 
radius x(t turn ) = v k i ck (to/ti) 2 ^ 3 tid. Where d = 0.8638 

As we shall see from the numerical results, the solution rapidly enters 
a self-similar growth phase, with turnaround radius scaling as x turn ~ t 2//3 . 
However, there is an intermediate time when the wake 'binds'. During this 
time, the second turnaround radius decreases, then increases once more at 
the onset of £ 2 / 3 growth. We can calculate the radius in the growing phase 
(after binding) using the approximation outlined above, but now, we consider 
a particle falling in from one side of the wake, crossing the symmetry axis, 



and flowing out the other side of the wake. This calculation is an extension 
of the Zel'dovich approximation into the non-linear regime since the particle 
trajectory we consider will cross other particle trajectories. However, the net 
effect of the gravitational field should be the same, as long as the particle 
turns around outside of the rest of the particles in the wake. Keeping only 
growing mode terms, we find the turnaround radius 

2 2 

Xturn = -fVkicklj-j U(-^-j (12) 

Including damping mode terms shifts /. Simulation results given below give 
a good match to this expression with / = 6.7. 



3 Numerical Scenario 

We use standard particle-in-cell (PIC) techniques to simulate dark matter in 
an expanding universe. PIC techniques use particles which statistically rep- 
resent mass in a fluid, in our case the fluid is collisionless. The particle masses 
are interpolated onto a mesh giving a mass density. Next, the potential is 
solved for by first Fourier transforming the density field, then solving an al- 
gebraic equation for the potential in Fourier space, then retransforming back 
to configuration space. The discrete derivative of the potential then gives 
the gravitational force, which is interpolated back to the particle positions. 
Once the force is known on each particle, a staggered leapfrog algorithm up- 
dates the particle locations and velocities and we continue updating particle 
positions and velocities in time. We check that energy is conserved by our 
code to a fraction of a percent. 

As initial conditions for the simulation, we give an ingoing velocity to 
two matter streams on either side of a symmetry axis and we set periodic 
boundary conditions. Our simulation has typical grids of 2048 zones. For 
cold dark matter, we use the same number of particles as gridzones. However, 
for hot dark matter, in order to reduce shot noise from the thermal velocity 
distribution, we use 512 particles per zone. This reduces spurious statistical 
fluctuations in the overdensity to Sp/p ~ 0.05. 

For hot dark matter, the physical thermal velocity peaks at t eq and then 
decreases as the scale factor during the matter dominated era \=j. The 
mean physical thermal velocity at t eq is given by the neutrino mass, which is 
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determined by requiring that neutrinos make up the critical energy density 
for an Q = 1 universe. We find 

v eqp = T vm jm v ~ 0.05 (13) 

where v eQp is the physical thermal velocity at t eq , T vm is the neutrino 'tem- 
perature' at t eq and m u is the neutrino mass. Thus, in hot dark matter 
simulations, a Maxwell - Boltzmann distribution of velocities peaked about 
this thermal velocity is then added on top of the streaming velocities of the 
particles. 

Throughout the paper, we assume a critical universe with Q = 1 and 
h = 0.5. 



4 Simulation Results 

We have checked the infall velocity by comparing with the infall velocity 
from the Zel'dovich approximation. A plot of the ratio of the Zel'dovich 
approximation velocity with the numerical velocity is given in figure |l|. Notice 
that the maximum difference is by a factor of 0.4%. 

As mentioned above, the analytical result for the 'second turnaround' is 
exact until just after the point where a particle crosses the path of another 
particle. In figure |2], we plot the analytical result versus particle trajecto- 
ries. The discrepancy can be seen once the particle falls back toward the 
symmetry axis. It is interesting to note that the wake width does not grow 
as £ 2 / 3 until after 'binding' occurs. The solution follows a transient from the 
formation stage to the growing stage, resulting in a smaller wake than would 
be predicted if one assumed the onset of self-similarity directly after the first 
particle falls back toward the symmetry plane. We also plot the growth 
(~ a{t)) given in equation |l!] from the point where this second stage in the 
evolution of the wake takes over. In figure |^ we plot the second turnaround 
radii (analytical and numerical) for the evolution of the wake until late times. 
The analytical curve matches well with the numerical wake radius relatively 
quickly (notice that the plot axes are logarithmic). If we extrapolate the 
curve to time t ? we find a wake size of about 2x turn = 0.55/i^q 1 Mpc for 
Zi = 10000 and v s = 0.5. 

In figures |] and |^, we plot the phase space of the particles in a cold 
dark matter simulation with Zi = 10000 and vi = 0.5c. The velocity axis 
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Ratio of Zel'dovich Approximation to Numerical Streaming Velocity 
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Figure 1: A plot of the ratio of the analytical infall velocity to the numerical 
result. After an initial correction, the infall velocity is correct to %0.4. 
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Particle Trajectories in Cosmic String Wake 
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Figure 2: A plot of the trajectories of a number of particles compared with 
analytical results. In the first stage of growth of the wake the leading parti- 
cle's trajectory is enveloped by the analytical trajectory for a particle being 
kicked away from the symmetry axis. In the second stage, where particle 
turnaround grows as a(t), the turnaround radius is enveloped by the analyti- 
cal i 2 / 3 power law. In this figure, solid curves with + signs indicate analytical 
results, other solid curves are numerical results. 



particle turn-around radius 



-1 



-2 



w 
o 
a. 






£ -6 



c 
o 

w 
o 
a 

o 



-9 



-10 



-11 

-2 



I I I I 


\ N 


^b, -t 


^^ + 


Tfcl + 


I I I I 



10 



log(redshift) 



Figure 3: Here, particle turnaround (ie. the width of the wake) (dash-dot 
line with plus signs) is plotted versus the analytically calculated growth law 
(solid line with circles and crosses) well into the evolution of the wake. Note 
the close correspondance over the final 3 orders of magnitude in redshift. 
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scales as a -1 / 2 and the position axis scales as a(t) in order to make the self- 
similarity evident. The phase space is similar to that shown in Fillmore and 
Goldreich JT5 in their figure 2. The main difference comes from the velocities 



outside the non-linear region. In the phase space of a cosmic string wake, the 
streaming velocity outside the non-linear region is constant. This is because 
the velocity perturbation from the string is due to the conical deficit angle 
which extends approximately to the Hubble radius. Thus, in the Newtonian 
approximation, all particles have the same velocity, regardless of distance 
from the symmetry axis. 

In figure |6], we plot the phase space of the particles in a hot dark matter 
simulation with Z{ = 10000 and v% = 0.5c. The phase space looks similar 
(except for a certain raggedness, which is from shot noise) to that of the 
cold dark matter wake at late times. In simulations, we find that accretion 
does not begin until the free-streaming length is smaller than the width of 
the wake. Before this time, the effective pressure of the neutrino gas keeps 
structure from forming on small scales. 

The neutrino free-streaming length grows until t eq , then decays as a{t)~ l l 2 . 
We can approximately calculate the time when wakes formed at t eq should 
begin to accrete by finding the time when the particle turnaround (calcu- 
lated above for cold dark matter) is equal to the free-streaming length. This 
calculation gives a result 

t accrete ^4.5xlO^-^—( b M. (14) 

bv kick \ t J 

Thus, accretion should begin at z accret e ~ 30. 

Numerically, we define a wake to be accreting once the overdensity Sp/p = 
1. We find hot wakes formed at t eq beginning to accrete in our numerical 
simulations at z accret e — 30 and second turnaround achieved at z turn ~ 8. 
After second turnaround the phase space of hot wakes resembles that of cold 
wakes with the overdensity rapidly becoming self-similar and growing roughly 
as a(t). We cannot make more quantitative statements due to the growth of 
spurious perturbations from statistical fluctuations in the initial conditions. 
However, after the onset of accretion in hot dark matter wakes, the evolution 
proceeds much as it would for cold dark matter with a rapid approach to 
self-similarity 
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phase diagram of cold dark matter wake 
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Figure 4: A plot at a late time (z = 7.202) of the particle phase space of 
a cold dark matter wake formed at Zi = 10000 with Vkick = 0.5c. Each 
turnaround in the phase space trajectory corresponds to a peak in the wake 
overdensity. 
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phase diagram of cold dark matter wake 
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Figure 5: A plot at a late time [z = 0.922) of the particle phase space of a 
cold dark matter wake formed at Zi = 10000 with Vkick = 0.5c. Notice the 
self-similarity with the previous figure. 
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Phase Diagram of Hot Dark Matter Wake 
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Figure 6: A plot at a late time (z = 7.212) of the particle phase space of a hot 
dark matter wake formed at Zi = 10000 with Vkick = 0.5c. Each turnaround 
in the phase space trajectory corresponds to a peak in the wake overdensity. 
The position of only one in 512 particles is plotted. 
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5 Discussion and Conclusions 

In this work, we have presented results from a high-resolution numerical 
simulation of the evolution of the hot and cold dark matter distribution in 
cosmic string wakes. We have also extended the Zel'dovich approximation 
and found expressions for the size of the dark matter overdensity in wakes. 
We have compared these and previous results derived using the Zel'dovich 
approximation to our numerical results and find good agreement. 

The most massive wakes today are those that have the longest time to 
accrete. These wakes are formed at t eq . Our simulations give a wake width 
of 1.17 h^Mpc for wakes formed at z eq . Hot dark matter wakes are kept 
from accreting until the mean free path of the neutrinos falls below the size 
scale of the wake. The redshift that this occurs for a wake formed at t eq with 
Vkick = 0.5c is z accre te ~ 30. We find good agreement with this prediction 
in our hot dark matter simulations, with accretion beginning at the time 
calculated from rough arguments. 

In a cold dark matter wake, our results show two distinct stages of be- 
havior of the wake width: the first stage is a brief epoch of essentially linear 
growth and contraction immediately following wake formation, and a second 
stage in which the overdensity begins to grow as particles accrete onto the 
wake. The transition between the linear stage and the growing stage is an in- 
teresting non-linear phenomenon. The fast approach to non-linearity makes 
the cosmic string wake a useful testbed for n-body codes. 

To get an understanding of the cosmological implications of these results, 
we use data from network simulations. Cosmic string velocities from simula- 
tions are distributed with a broad peak in the distribution at about v s = 0.3 
on horizon scales. Thus, the average cold string wake formed at time t eq will 
form a wake of mass 3.2 x 10~ 8 of the mass in the horizon at to- Combined 
with network simulation results (mentioned above) that there are ~ 10 — 30 
strings per horizon volume, all of the matter in a cold dark matter universe 
would be non-linearly accreted at present. Some caveats to this result are 
that our simulations assume a matter dominated universe, and that coarse- 
graining in network simulations may lead to an incorrect estimate of the 



number of strings per horizon [ISfl. Nevertheless, a significant fraction ~ 1/6 
of the matter in the universe would still be non-linear even if there were only 
~ 1 string per horizon. 

In the analogous hot dark matter universe, significantly less matter will 
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be non-linear due to the late accretion time for hot wakes. Hot wakes formed 
at t eq by cosmic strings with velocities v ~ 0.5c begin to accrete at z ~ 30. 
Early structure formation is therefore not a problem in the HDM plus cosmic 
string model (strings with higher velocities will form wakes that accrete even 
earlier) . If we use the cold dark matter expressions for wakes forming at this 
redshift, we find that a fraction 0.002 of the matter in the universe is in hot 
wakes. This fraction will go up if we integrate over all wakes. Also, recent 
work on the string network scaling solution shows that the length in long 
strings is roughly two orders of magnitude higher in the radiation era than in 
the matter era li-il. This work shows that there is a relatively long transition 
between the two scaling regimes, and thus, there will be more long strings 
than previously thought around t eq . This effect will again push the amount 
of accreted matter upwards. 
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